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ABSTRACT 

Performance of some fibre optic digital communication 
links with different source- detector combinations and fibre 
lengths are evaluated. Fibres are dispersive in nature and 
spread the optical signal on transmission through them. 

This dispersion with additive noise cause intersymbol inter- 
ference (ISI) at high data rates, A decision feedback 
equaliser (DFE) has boon used to reduce the effect of ISI* 

To evaluate the performance of the DFE, the optical 
systems with different lengths and sources have been 
simulated on the digital computer DEC- 1090 system. Simula- 
tions have also been carried out at much reduced data rates 
without the DFE. Results indicate that similar performance 
in terms of the probability of error of the system can be 
achieved without an equaliser, but, at much reduced data 
rates. P(e) bettor than 1.25x10'" 7 has boon obtained with 
laser source and p(e) of 2. 5x10" 5 with LED source for a 
SNR of 14 dB. 
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CHAPTER 1 


INTRO DUCTIOi'I 

1.1 GENERAL: 

The use of light for senciin> : message information 
goes far back into time. In 1876 Alexander Graham Dell [l] 
explored the use of light for nucHo Communication. The 
invention of the LAoER in 196U--61 stimulated widespread 
interest in the use of light for communication purpose. 

In 1966 K.C. Kao and G.A. Ilockhein [3] proposed that optical 
fibre be used as dielectric wave-guides in under ground cable 
to carry information. The advent of optical fibres added 
a new dimension to optical communication. In I960 typical 
fibre losses were above 1000 dB/km, but by 1970 losses 
under 20 db/km wore achieved. 

Presently installed fibres have losses in the range 
of 3 dB/krn or less at wave lengths of 820-850 nm. Laboratory 

) 0 V/ 8 $ 

samples have shown losses as/Ci.d? dB/km at 1550 nm wavelength 
The present objective is to develop fibres with losses around 
1 dB/km for commercial use. ! 

The advantages offered by fibre optical com,! uni cat ion 
are enumerated below. 
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a) The system is immune to In. J , RFI and cross 
ta Ik , 

b) It offers large bandwidth and low loss, 

c) Small size and weight with consequent ease of 
installation and reduced transportation volume. 

d) Potential low cost 

e) No electrical ground loop or short circuit problem 

f) Reduction of total system povv'er consumption. Actual 
power output used Is of the order of mw. 

1.2 OPTICAL FIBRE COUUNICATION SYSTEM: 

The block schematic of the fibre optic digital 
link is shown in Fig. 1.1. The system can be divided into 
two main parts - one part which deals with the electrical 
signals and. the other part dealing with the optical energy. 
The later part is shown within the dotted lines in the 
figure and includes the optical source, the optical fibre 
and the photodetector. 

The different subsystems in the fibre optical digital 
link are briefly explained below. 

Dri ver: Like in any communication system the input signal 
or data first modulates a carrier which is propagated 
through the channel. Here too the input modulates the optica 
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source directly, when the input .is 1 1* , the optical source 
is turned on otherwise it is hept off. 

Ele ctr o- optica 1 Transduce r : oOi.n concluc tor laser and light 
emission diode (Li-D ! s) are being" used for converting 
electrical signals into optical signals. In the system 
under discussion we have used both, the laser diode as 
well LED. 

.Fibrje: The optical signals are propagated through 

the optical fibres channel. Those signals undergo some 
losses and distortion depending upon tho characteristic of 
the fibre. The distortion is such that it spreads the 
signal pulse in time domain. This is mainly caused by two 
types of dispersion in optical fibres, i.e. intermoda 1 and 
material dispersion. However, in graded- index fibres 
only tho m a ter i; 1 1 d is p o r s i on i s s 1 gn i f i ca n t . 

PjL^iPA { t. c .^ r A c - The optical signal is received 

by a pho to-U'O tector , which converts the incident optical 
signal into an electrical signal. Tho photo do tec tors 
mainly of two types. 

1. P-1--H Diodes 

2. Avalanche photodiodes (,PD's) 

Fro nt End Receiver : The front end receiver circuit is a 
combination of a pre-amplifier and a post-amplifier. As 
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the data rate is very high, the pr e. .amplif ier is implemented 
using a FET, CA-3127E. 

Equalizer ; The distortion in the optical fibre channel 
spreads the signal pulse in time domain. The material 
dispersion and the additive noise of the system cause ISI. 

An equalizer is used to reduce the effect of the ISI. In 
the system under discussion we have used a Decision feed- 
back equalizer, which reduces the magnitudes of the inter- 
fering samples. 

Bit Synchroniser: The bit synchroniser recovers the symbol 
timing clock from the received signal and selects the 
sampling instants. The equalized signal is sampled at the 
right instants and at the rate at which it is received. A 
decision is then made about the campled signal using a 
threshold detector, and we get the required signal after 
it is decoded suitably. 

In any communication system the input signal modu- 
lates a carrier which is propagated through a channel. Here 
we consider the input as a Binary data sequence of 1 1' or 
*0* . The optical source output intensity is modulated by 
the input electrical signal. Thus when a * l * is trans- 
mitted the optical signal having some intensity propagates 
through the fibre. The impulse response of the fibre 
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spreads the signal pulse in time domain. The attenuated 
and dispersed optical signal when arrives at the receiving 
end of the fibre, it is converted into an electrical signal 
by direct detection in the photo detector. During the detection 
process certain noise inherent to the detection mechanism and tl 
front-end circuits got added 'to the* signal. The effect of the 
noise- depends upon the typo of the- photo detector used. This 
will be discussed in a la tor section. 

By the principle of superposition for optical 
fibre transmission (Personick [ L f ] ) single pulse descrip- 
tion can bo extended to obtain a model for the transmission 
of an entire sequence of ON and OFF pulses. The photo 
detector output r(t), can be written as 

r(t) = E g K s( t-t^) + n(t) (l.l) 

K 

whore the time instants t^ form a poisson process having 
intensity /\(t), with 

A(t) = E a jL p F (t-iT) + /\ Q (1.2) 

where c^. = avalanche gain factor 

s(t) = output pulse of photodetoctor 
n(t) = additive noise 

T = signalling interval 



7 


/ \ Q >_ 0 is the dark current 
a^ = 0,1 equiprobable data symbols 
Pp(t) > 0 pulse response of the fibre 

The output from the photodetector as given by 
eq. (l.l) is then filtered to reduce the effect of noise 
by passing through a low noise pro-amplifier. A decision 
is then made about the received signal using a threshold 
detector and wo get the required signal. 

An optical fibre communication system might differ 
from the one discussed in this section depending upon the 
source, data rates and fibre length used. These aspects 
have been discussed in the chapters to follow. 

1.3 NOISE IN PHOTO- DETECTORS: 

The photo- do tec tor current consists of the sum of 
displacement currents of individual hole and electron 
pairs generated within the photo-detector. The time at 
which electrons are generated duo to absorption of the 
impinging photons are random and obey Poisson distribution. 
Every individual optical pulse (correspondence to ’l') 
differs from the average by some unpredictable amount 
and this difference can bo termed ns a signal dependent 
noise. In practice electrons are also omitted when no 
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optical signal (photons) fall on the photo-detector. This 
is called the dark current and is modelled by introducing 
a constant additive intensity function before the 

detector (see eq. (1.2)). To transmit a 'O', the optical 
pulse is not propagated for that bit interval, but, the 
detector emits electrons due to thermal ionization (dark 
current) . 

Another source of noise is thermal noise which is 
always present in electronic circuits. A circuit model 
(Fig. 1.2) is shown to illustrate the various noise sources. 
The output current from the photo-detector is given by 



Fo <■ S'SO cs 2 A% 

COi'f 1 '- dCtJ 


fig) ±‘Z ^ETfcCTOR 0 4) £ L 




9 


i( t) = i°(t) + i d ( t) + i th (t) (1.3) 

where, N 

i S ( t) = 2 g i h ( t~z . ) 

j=l J S J 

N d (t,t+T) 

I d (t) = 2 g. h d ( t-z. ) 

j=l 3 a J 

c* 

whore i°(t) is the shot noise process (signal dependent) 
and N (t,t+T), the number of photo-electrons generated in 
a closed interval [t,tfT], is random and has Poisson 
distribution with moan /\ e (t) 

o 

t+T 

whore A ( t) = / p(t) dt 

s t 

z. is the time when the j th electron is emitted* h c (t) 
is the time invariant filter modelling the component 
function of the shot noise. i d (t) is the dark current 
process with an intensity / \ Q , 

l ( t) is the additive thermal noise process and 
gj is the gain of the photo diode. 

In the case of a PIN diode detector g.=l, and pulse 
signal-to-noiso ratio, /h s 2 (t)/N Q , is of the order of 
-20 dB. Therefore, the shot noise is neglected In compa- 
rison with the thermal noise in the case of PIN diode. 
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In the case of APD detector the gain is large. This 
increased g a , increases the shot noise. Pulse signal- 
to-noise ratio / g“^h ( to( t)dt/N_ is of the order of 20 dB. 
Here the shot noise is quite large compared to the 
thermal noise and hence thermal noise is neglected when 
APD detector are used. As we have used a PIN diode 
detector in the F.O.C. system studied, we can neglect the 
shot noise component, 

1.4 EQUALIZATION; 

The model of the F.O.C. system has boon discussed 
in Section 1.2, In that section we have noticed that 
received pulses are spread out in time. Therefore, the 
transmission of higher and higher data rates will be 
possible if the fibre optic medium accomodates larger BW. 
Otherwise, with the higher data rates, ISI will increase 
and in turn increasing the probability of error. 

There are two general methods for equalization; 

a) Frequency domain equalization, and 

b) Time domain equalization 

Here wo will discuss Time-domain equalization. 
Since early i960, there has boon a lot of interest in 
time domain equalizers. 
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A typical fibre optic receiver for digital signalling 
schemes consists of a photo- do tec tor , an amplifier, an 
equalizer and a decision device, along with the necessary 
synchronization circuitsas shown earlier. Vvc treat the 
fibre optic channel as linear in power [4,5], i.o. 

CO 

P(t) = £ X ( i ) h ( t-iT ) (1.6) 

i=-co iJ 

whore , 

P(t) - is optical power driving the photo- detector 
I(i) - optical power from source 
h ( t) - is the fibre response 

Jr 

The average signal at the equalizer output is 

CO 

< y(t)>= z I ( i) p(t-iT) (1.7) 

i=— oo 

where p( t) is the overall pulse response. The equalizer 
output is given by 

oo 

y(t) = E I(i) p( t~iT)+n( t) (1.8) 

i=_oo 

whore n(t) is noise, modelled as a Gauss ian-plus filtored- 
Poisson process. 

In time domain, the received signal after the 
photo-detector is fed into the equalizer which is a 
delay line with multiple taps, viz. Transversal filter of 
Fig. 1.3. The tap weight vector 0 ~ [C^C^. . • c n ]'*' is 
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fi t> 1*3 tappeo oetAV-L 




i n C a i i ^ C O 


computed in order to satisfy a suitable criteria. Two 
equalization criteria are often used, viz. zero-forcing 
as suggested by Lucky [6] and minimization of the mean- 
square error (MSE) [7], In the zero forcing criterion the 
output distortion ’ D 1 is minimized, whore 


D = 





- v i, 
K=~oo, 0 


and the peak distortion of the unequalized pulse < 1. 

In first case, C is computed by solving a linear 
system. The second criterion minimizes the MSE between 
the received sequence estimate and the transmitted sequence. 
From a mathematical view point, this approach searches for 
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linear least squares solution. The main draw back of the 
tapped delay line equalizer results from the fact that the 
length of the equalizer tends to infinity. 

In 1967 Austin [8] determined a new suboptima 1 
receiver structure which is non-linear and uses its previous 
decisions to reduce ISI in high speed data transmission. 

It is called a Decision Feedback Equalizer (DFE) as it feeds 
back the output symbols through another TDL for comparison 
with the forward TDL output. After this a lot of work has 
been done on DFE. The main advantages of DFE are: 

i) Number of taps required for the forward TDL is 
much less than the transversal filter equalizer 
(unitary equalizer) tap delay lino. 

ii) The number of taps required for the backward TDL 
is also limited. 

iii) The performance of the DFE is better than the 
unitary equalizer for the same data rate. 

1.5 EQUALISATION IN FIBRE OPTIC DIGITAL TRANSMISSION: 

Since the development of fibre optic digital commu- 
nication, a lot of work has boon done for the minimization 
of ISI and noise. 

Foschini ot.al. [9] in October 1975 worked on 
optimum direct detection of digital data signal. In their 
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approach they processed the- output of the photo-detector, 
so that the- probability of error is minimum using the 
Vitorbi algorithm detection principle. Vitorbi algo- 
rithm detector is difficult to implement in hardware for 
high data rate systems with present day technology. A 
sub-optimal solution using a DFE as the first stage in 
place of the first VA detector of above has been considered 
by Govind Sharma [10]. The second stage of the receiver 
as obtained by him is also quite complex as it requires 
large memory and a very high sampling rate. 

Doqliotti et.ai. [ll] in 1976 considered the 
problem of base-band equalization in optical PCM fibre 
system. For the given data and fibre characteristics they 
have suggested an equalizer transversal filter, where 
the tap weight vector C % [C^ , ... C fi ] are computed in 
order to satisfy a suitable criterion for the minimization 
of ISI. In this paper as well as one by Dogliotti and 
A. Luvison [12], the authors have minimized IS 1 by sequence 
estimation technique. In such a case, the received process 
is written according to its 'state equations' and the 
optimum estimator is derived as a Kalman Filter [13], For 
TDL structures, the MIvlSE criterion guarantees bettor 

performance as shown in Fig. 1.4. This figure corresponds 

_g 

to a probability of error of 10 and shows that A/MSE 
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equalization is better than zero forcing equalization and 
also when o/T<l, then we can even ignore the equalizer 
provided that SNR is high enough, 

D. G, Messers chmitt [14] considered a digital fibre 
optic transmission system with Poisson signal statistics 
and additive wide sense stationary noise. The criterion 
of optimization is MSE between the decision threshold 
input and the current data digit EL. He has considered 
the zero forcing criterion as well as MMSE for both linear 
and non-linear (DFE) equalizer and obtained solution for 
each criterion. Some special cases like white noise and 
band limited noise have also been considered. Capt. R. 

Wahi [15] designed and evaluated the performance of the 
DFE based on the Austin [8] sub-optimal receiver structures. 
Ho also included predistortion filter in the transmitter and 
then evaluated the performances of system. At lower SNR 
values pre-distortion filter reduces error and gives a 
gain of approximately 1.5 dB. G, 'famburclli [16] has 
examined the performances of decision feedback and feed- 
forward receiver whoiu the post-cursor correction is 
carried out only once while the precursor correction takes 
place many times without requiring complex circuits as 
the relevant transversal filter consists of only one or 
two taps. He has shown that the post-cursor of the 
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incoming signal is compensated by the decision feedback, 
thus at the output of a first decision circuit the trans- 
mitted signals arc evaluated with an error probability 
relevant to a decision feedback receiver, and precursor 
are compensated by other decision circuits which follow the 
first. He lias shown that this arrangement is superior to 
the existing ones and allows a transmission rate equal to 
twice the Nyquist rate. 

1.6 PROBLEM DESCRIPTION: 

The problem considered in tnis thesis is to evaluate 
the performance of the Digital Fibre Optical Communication 
System at different data rates. The following systems are 
cons id or ed. 

a) Source-laser diode, link length 2.155 km 

b) Source-LED, link length 2.155 km 

c) source-laser diode, link length 4.310 km 

At higher data rates, noise and I Si affect the performance 
of the system. This can be reduced with the help of an 
equalizer. W e have used' a DFE for this purpose. The tap 
weight values for the DFE arc calculated by Monson' s 
Adoptive linear filtering method. 
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The F.O.C. system is simulated on the digital 
computer DEC— 1090 system to find the maximum data rate 
that can be transmitted over the system with and without 
the DFE for different SNR values, 

1.7 ORGANISATION OF THESIS: 

This thesis comprises 5 chapters. In Chapter 2, we 
have briefly reviewed the work done on adaptive equalizers 
by Peter A/tons on. Chapter 3 describes a fibre optic commu- 
nication system and its link parameters and the design of 
the DFE. It also describes the impulse responses for 
various link lengths and the bandwidth of the channel. 

Chapter 4 describes the simulation of the DFE on the 
DEC 1090 system. Tap gain settings and the probability of 
error for various SNR’s at different data rates are tabulated. 

Chapter 5 concludes this thesis with a discussion of 
the results obtained, and some suggestions for possible 


future work 


CHAPTER 2 


DECISION FEEDBACK EQUALIZER 

2.1 INTRODUCTION: 

Synchronous, digital communication systems are 
degraded primarily by noise and intersymbol interference. 

To combat these effects , equalization techniques [l7],[l8], 
have been suggested. In 1967 Austin [ 8] considered a 
feedback receiver that compensated for the effects of inter- 
symbol interference due to both future and past symbols. 

His optimum structure consists of a ma tched-filter , transversal 
f il ter combination (forward filter) operating on the 
received signal and a feedback linear filter operating on 
the past decisions. The past decisions have been assumed 
to be correct. Non-s tationary effects of the channel were 
not considered. 

Peter Monsen [19] has worked on adaptive equalization 
again under the fundamental assumption of no decision error 
so that the output digits can be assumed correct. The i 

optimality criterion is minimum-mean- square error (MMSE), ; 

Under this criterion, the receiver is shown to be superior 
to one without the feedback path. ; 

This feedback receiver, besides mitigating noise and ISI 
acts to eliminate timing jitter and Doppler shifts in the ; 

, i 

channel. i 
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Some theoretical calculations and experimental results 
indicate that the error propagation effect due to incorrect 
decision is not large. (Austin has also shown it). Since the 
I Si span is normally short, a discrete Markov process model 
for analysis of the error propagation effect can be used to 
determine the error propagation effect. 

2.2 SYSTEM DISCRIP TION s 

The communication model under study considers trans- 
mission of serial binary data over a linear channel for which 
the dispersion has width of the order of the data rate 
reciprocal. The channel dispersion produces ISI which can 
be excessive. 


The communication model is shown in Fig. 2.1. If the 



A .1 F.OMi.i 2tk -Y/S tCI'u 
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channel Impulse response is h(t), the received signal r( t) 
is given by 


r(t) = Z S . h ( t~ iT) + n(t) (2.1) 

i 

where T is the interval between symbols. 

The binary digits 3^ are equally likely, and have 
a correlation function. 




( 2 . 2 ) 


The noise term is a realization of a zero-mean white Gaussian 
process with spectral density p . 


The receiver consists of tv/o linear filters, the forward 
filter that operates on the received signal, and the back- 
word filter that operates on the .reconstructed data. Delays 
in the channel impulse response can force the use of a 
noncasual forward filter when the transmission delay is nor- 
malized to zero. Accordingly, the impulse response of the 
forward filter y(-t) is defined over the entire range 

_ co < t < co. 

The backward filter has discrete Input and output 
and thus can be realized as a transversal filter with symbol 
interval tap spacings. Because the filter is located in a 

y * " I* 1 ' 1 * 1 1 ' mmmmmrnmmmmmmmm 

It is notations liy convenient to define the forward filter 
Impulse response as the function t(-t)* 
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feedback loop, it can only operate on past reconstructed 
digits* Its impulse response can then be written in terms 
of the transversal filter tap gain b^ as 

CO 

b(t) = E b. ' ( t~iT) (2.3) 

i=l 1 ' 

The reconstructed digits are formed by a decision device 
that produces the sign of the analog sample at the summed 
output of the two filters. 

2.3 EQUALIZATION FOR KNOWN CHANNEL: 

To solve the time varying digital communication 
problem, various optimization criteria will be investigated 
for a stationary environment. After considering the 
different criteria, adaptive techniques will be considered. 
The input of the forward corrective filter is the sum of 
the signal and noise terms. 

CO 

r(t) = E SA h(t-KT) + n(t) (2.4) 

K=~oo 1 

The forward filter linearly operates on the recon- 
structed discrete synchronous output. The outputs of the 
forward and backward filters are combined. Tho impulse 
response of tho backward filter can be written as 

OO 

b(t) = E fcc f; (t-KT) 

K=~<» 


(2.5) 
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Thus the combined signal into the decision device at the 
sampling time j T is given by 

' s. = Y. - z. = E y(u) r( jT+u ) du - E b, S- 

3 3 3 j=~oo (=1 3 " l 

( 2 . 6 ) 

/\ 

Under the assumption = S^, the input to the decision 
device can be split into the signal, noise and intersymbol 
interference terms. These components found from equations 
(2.4) and (2.6) are as follows 


S j' signal = S j X' - h -o 

U I intersymbol ~ U t X'~l 
interference 

I noise = X' G.j 


(2.7) 


E S . , 

=1 ' 3 • 

( 2 . 8 ) 

(2.9) 


A receiver design depends upon the criterion which 
minimises the intersymbol interfor mco term and the noise 
term, 


2,3,1 Unitary Equalizer; 

In this section following types of unitary equalizers 
are described in brief. 


1. Transversal Filter Equalizer 

2. Zero Forcing Equalizer 
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2. 3. 1.1 Transversal Filter Equalizer: 

For a unitary equalizer in a receiver consisting of 
only a forward filter and a decision device, the inter- 
symbol .interference term can be eliminated if a function 
Y ( t) can be found, y( t) may be a tandem configuration. 

OO 

Y(-t) = / g(-u) \J( t-u) du (2.10) 

— CO 

where g(~t) is the impulse response of a fixed filter with 
bandwidth roughly equal to that of the channel, and W( t) 
is the impulse response of a transversal filter with tap 
spacing equal to the symbol interval. 

If the transversal filter lias N negative - index tap 
outputs and ^ positive index tap outputs, W( t) takes the 
•f orm 

n 2 

h(t) = E W. f, ( t-iT) (2.11) 

i=-N x 1 ' 

where is the ith tap gain. These tap gains are then to 
be chosen to minimise the effects of intersymbol interference. 
The forward filter impulse response is given by 

N 2 

y(-t) = E Vi, g(iT-t) (2.12) 

1 =-^ 

and the system pulse response is 
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qv- 4 q(KT) = Y' ^ ^ V/. 

R i=~N 1 1 

Tho integral in eqn. (2.13) is 
input to tho transversal filter 


CO 

f g(u+ iT)h(u+KT)du 

«~«CO 

(2.13) 

the pulse response at the 
at the sampling time (K-i)T. 


Let 

oo 

9 V = f g(u) h(u+KT)du 
K 


(2.14) 


Thus eqn. (2.13) can bo written as 

N 2 

q K = * w i e K-i = S’ e K 


where, 


w = 


NO 

/ - 

il 

® K+N 1 


„ 

I'l- 1 




2 i 






(2.15) 


( 2 . 16 ) 


For tho unitary equalizer tho intorsymbol interference term, 
eqn. (2.8), is then 


3 


inter symbol I ^o 
in terfer once 


(2.17) 


with a finite number of degrees of equalization freedom. It 
may not be possible to force eqn. (2.17) to be identically 
zero. A possible optimization problem is 
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min £ cP: , q = constant (2.18) 

> / K 0 

* £ 0 

The performance of the unitary equalizer can be 
consideroci by first evaluating the performance when and 
N 2 are infinity and then investigating the case when or 
Np are restricted to a finite number of degrees of freedom. 
V/hen and 3x0 infinity, the ISI can in most cases be 
reduced to zero if 

OO 

% = f W i S K-i = °> K * 0 

X— — 00 

2. 3.1.2 Zuro Forcing Equalizer: 

Intersymbol interference can also minimised by a 
zero forcing equalizer as suggested by Lucky [6], The 
peak distortion for the pulse q(t) is 

co 

D = to kSL. |c,kI (2 ' 19) 

The input pulse to the equalizer will bo assumed to have 
initial peak distortion D . The reference sample x Q win 
also bo normalized to unity. The problem is to determine 
the tap gain coefficients for a tap transversal equalizer 
which minimises the final peak distortion. 

Lucky establishes that at minimum distortion, at 
least L~1 values of q^ must bo zero. This gives L-l equations 

N 
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in L-l unknowns, but the question remains as to which L-l 
zeros in the output response does one force to achieve 
minimum distortion. A key theorem is then proved which 
shows that if the initial distortion 

D 0 = -ra-T s I e 0 l (2 - 2 °) 

* 0* K— OO 

K^o 

is less than unity then the minimum output distortion must 
occur for those L-l weight values which simultaneously 
cause the corresponding output response sample to be zero, i. 

q K = 0 *(£ K i 0 (2.21) 

when D q > 1, the settings which force corresponding output 
zeros may not be the optimum values, and a distortion greater 
than the minimum may result. 

This approach is designed for channels where not only 
is intersymbol interference the limiting factor (since noise 
power is not considered) but, also whore the intersyrribol 
interference is bounded above. This situation exists to a ; 
large extend on telephones channels where the technique 
has boon applied, 

2.3.2 Feedback Equalizer: 

Last section described the class of unitary equalizers 
which operate on only the received signal. The feedback ! 
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equalizer also considers the reconstructed data sequence 
o^, which has been assumed to bo the transmitted data. This 
section describes the choice for forward and backward 
filters for an intersymbol interference minimize Lion , The 
prefilter is taken as a matched filter so the input noise 
power to the transversal filter is minimized. 

From eq. (2.C) the ISI term at the input to the 
decision device is 



= 2 

inter symbol f^o 
in ter ferenc e 


CD 




( 2 . 22 ) 


where ^ 

q= ? (2 ' 23) 

l 

For minimization of in ter symbol interference, the feed- 
back filter bo 


1 + 

O v 

i'v 


AS. 


w, 0, 1 , 

3 i^“j 


is. 


JL f ^ f * m • f 


( 2 . 24 ) 


For this choice of a feedback filter, all the ISl due to 
past pulses is concolcd exactly. The ISl power becomes 


>) \ "1 
Y- E 

1 -=“°° 


% 


a 2 


j*-K 


3 


( 2 . 25 ) 


which, with an infinite number of degrees of freedom can be 
reduced to zero. Tho optimum weights solve the semi- 
infinite set of equations. 
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CO 

_E W t + 0 K ^ ± = 0, K < 0 (2.26) 

X=;«»oo 

Peter jvionsen [i9 ] , solved these Gelations by spectrum 
factor iza tion. 


Since unitary and feedback equalizers reduce the I SI 
to zero, one would quaj.ita ti vely expect the feedback 
equalizer to be superior because of its noise advantage. 


2.3.3 Mean- square Error minimization : 


For channels which arc not limited solely by 
either I SI or additive noise, an optimizing criterion should 
be chosen which minimises both the effects. The feedback 
receiver minimises the mean square error. The error is 
defined as the difference betv^een the analog voltage at the 
input to the decision device and the corresponding trans- 
mitted binary digit, i.e. 



(2.27) 


where 




t=l 






bj 


Y* n^ 


(2.28) 


Here error is produced by both 1 31 and additive noise. The 
minimum rnean-square-error equalizer minimizes the mean 
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square of eqn. (2.27) with respect to both forward and 
backward filters. By including the backward filter in 
the optimization. Monsen has shown that the jVJviSE feedback 
equalizer is superior in the absence of decision error 
to the j'viiwSE unitary equalizer. 


This criterion of min iraun>- :oan-squarG error (MMSE) 
has been compared with the minimum probability of error 
(mPE) criterion. The later is the boot criterion for the 
class of receivers that consider only one symbol at a time. 
The APE receiver and its associated probability of error 
were determined by using steepest descent techniques. It 
was seen that the j\p£ receiver h..d siqnif icantly better 
performance than the J-jviSE receiver. However, it was 
finally concluded from this study that although the two 
criteria are not completely equivalent in performance, a 
vast majority of situations do load to equivalent perfor- 
mance. The simplicity in the use of the MMSE criterion is 
then the deciding factor in the choice of the criterion. 

The minimisation problem to be studied can be 
stated as 


min iviSE = 

X, b i 


'If ” n 




■ - f ■ 


(2.29) 


The input to the forward filter In the vector notation will 
be defined as = r(tfMT), so that 
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t 


K 


- . r k X. “ J ^ b t 


Then let 


(2.30) 

(2.31) 


HK 



1 — 1 
i 

f 

h 


1<-2 


b 2 

== 

* 

, b = 

j 

o 


0 

; i 


• 


(2.32) 


be semi- in finite column voctorc, co the decision device 
input is 


% = - r i x - ft 


ct 


b 


(2.33) 


The jv'jviSE optimization is invariant to the choice of time 
origin so that for simplicity the zeroth sampling instant 
is chosen. The minimization c^n now be cast into a more 
fa mil ia r f or m , i . e . , 


min j'.GE = ( r ’ y_ 3’ b -3)" (2.34) 

X>ll 

By the introduction of partitioned vectors, oqn. (2.34) can 
be cast in the functional form of a single vector variable, 
as follows; 


min MSB = (f« V. 
f ‘ ' ' 1 


,S K 


) 


2 


(2.3b) 
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where 



The solution of ecjn, 
set of- equations. 


V K 


-K 



(2.36) 


(2.33) leads to the solution of the 


< v,,. V’ > f = <V K 3,, > (2.37) 

K an. K 

where the < > brackets are employed to denoted the 
expected value. 

The solution for the forward filter impulse response 
is found to be of the form 


Y(--t) = E W.h(iT-t) (2.38) 

i 

which 5s the representation for a matched filter in 
tandem wi th a transversal filter with symbol- interval tap 
spacing. The transversal- filter impulse response is 
written in terms of the tap gains as 


W(t) 


The receiver is 
as socia ted wi th 
a sso cia ted with 


5 w iS (t ~ iT) 


(2.39) 


then specified by the tap gain set^WjJ i=“‘ 1 to 


the forward filter and the tap gain set i b^ji=l't 


i H 


the backward transversal filter. The 


solution for these tap gains were d-termined by monsen in 
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terms of their ' Z 1 transform. The technique of spectrum 
factorization £>5] was employed to determine the solution 
for W(z) and B( z) . The solutions are 

W(z) = [U”( z) U + (o)] -1 (2.40) 

B(z) = T[U + (z)/U + (o)' (2.41) 

The positive superscript indicates a spectrum factor with all 
its poles and zeros outside the unit circle, whereas the 
• negative superscript factor has them all inside. Since 
Vi(z) has no singularities out side the unit circle, it 
follows that the forward tap-gain function is anticausal 
i.e. 

W. =o i > o 

The forward transversal filter thus operates only 
on future values of the matched filter output. 

Thus, in the absence of decision errors, the inter- 
symbol interference dual to past digits is cancelled 
exactly by the backward filter if the source digits are 
uncorrelated. Under this condition the forward filter then 
minimizes the mean-square error duo to noise and the ISI 
due to the future digits. It is interesting to compare 
the optimum filter and MMSE for the feedback and non feed- 
back cases. The non-feedback or unitary equalizer (UE) is 
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a special case of the feedback equalizer (FE), namely, 
when 8(f) is set identically to zero. Thus the optimum 
FE is superior in a JviLSE sense, to the optimum UE. 

2.3.4 a Gradient Technique Realization: 

For the fixed known channel, a s teepest-descent or 
gradient technique was found for the Iol minimization 
equalizer. For the ivimSE receiver, a gradient technique 
would need to synthesize three separate filters, namely, 
the matched filter, forward filter and the backward filter. 

In any real system the channel will be band limited 
effectively, so the channel and the matched filter can be 
s pec if i ed in sa mpie d— da ta f oxm. 

The received signal r(t) in sampled data form is 

CO 

r(tfKT) = S r[ | + KT] sine (Ot~i) (2.53) 

The optimum solutions are identical to the one obtained for 
the mean Square error in Gcction 2.3.3. The difference 
is that the forward filter can now be synthesized by one 
transversal filter with tap s pacings equal to the Nyquist 
interval (!/•>). 


A gradient algor it Inn for a i.j»SE receiver takes 
the form 


f 

-n+1 


■n 


a [ V V 1 f 

. - --ft 


~n 



(2.54) 
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where f is the receiver vector (2.36) at the nth 
-in v 

iteration. Eqn. (2.34) can be c; ct into the form : 


r , = x 
••n+1 -n 


«V G n 


(2.55) 


where e is Lhe error at the nth iteration. The gradient 
technique and the adaptive receiver implementation using 
this techniques are discussed in the next chapter. 

2.4 ERROR PROPAGATION EFFECT: 

The DF£ uses receiver decision in an attempt to 
cancel 131 due to the past pulses, A fundamental assump- 
tion has been made that the reconstructed data sequence 
is identical with the transmitted data. In reality it is 
not true. Decision errors will occur and they win cause 
two effects. The first is that of an adaptation error. 

Since the reconstructed data sequence is employed as a 
test signal to determine the learning algorithm parameter. 
Decision error will further perturb the ideal equalizer 
tap gain setting. The second effect is more severe. 

The DFE decisions are linearly filtered in order 
to develop the appropriate voltage for the ISI cancellation. 
When the decision is wrong, the probability of decision 
error in the immediate future is enhanced as the decision 
error propagate through the backward filter. This effect 
is thus termed Error Propagation. The resulting error 
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propagation will be larger than the error propagation under 
the No decision error premise. 

Mon sen [20] analysed the magnification of the error 
probability due to the Error propagation effect for a 
single echo channel and a truncated geometric response 
and observed that degradation due to error propagation 
effect is very small. 

Since the error propagation effect is negligible, 
the DFE is superior to the linear equalizer. Moreover, 
because the backward filter eliminates interference due to 
all the past digits, the number of transversal filter taps 
in the forward filter can be reduced. This leads to an 
implementation advantage and improved adaptation as the 
self noise due to adaptation increases with the number 
of taps. 



CHAPTER 3 


D.F.E. FOR FIBRE OPTIC LINK 


3.1 INTRODUCTION: 

An important area of application for optical fibres 
is the transmission of high speed digital data signals. Even 
though the dispersion effects caused by a fibre are greatly 
reduced with respect to other kind of channels, the main 
limiting factor in increasing the data rate is imperfect 
channel characteristics. 

The approach to match the exact characteristics of 
the system to reduce the effect of noise and I si is not 
practical since it is difficult to design transmitting and 
receiving filters having the exactopt imal characteristics. 

Vie have seen that the DFE is a simpler alternative which 
provides a considerable reduction of the effects of noise 
and ISI. In the next section adaptive DFE based on the 
gradient technique has been described. 

3.2 ADAPTIVE DFE: 

For a fixed known channel, an equalizer realized 
by gradient technique minimises the mean-square error. 

This technique is used to synthesize three separate filters, 
i.e. the matched filter, the forward filter an-d the backward 
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filter. Furthermore, the matched filter impulse response 
is a continuous time function. In any real system the 
channel win be effectively band limited to some region 
| f j <_ 3/2. So the channel and the matched filter can be 
specified in sampled-dai. iorm. The continuous vector 

I / \1 OO 

X = i y ( t)f can then bo replaced by the infinite - 

L / t=— 00 

dimension dis crete- index vector. 


Y ± = Y(i/B) } “ 


(3.1) 


and j completely describes y( t) by the interpolation 
formula 

00 

y( t) = E Y n - sinc(Bt-i) (3.2) 

i=-oo 

If the channel is band limited to B/2, it is reaso- 
nable to position a filter prior to the receiver which 
will also band limit the additive noise to the range 
| f S < B/2. Then the received signal r(t) also has a 
sampled data representation of the form 


r(t+KT) = I r[~ + K.T] sinc(Bt-i) 

P 




^r(t*KT)i 


J i=— 


So the vector r, 

Ts. 

= ( r L ff + KT ]) 

The optimization problem is identical and if infinite-length 


CO 

Xrr-oo 


can be replaced by 
(3.3) 
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transversal filters arc used to realize the forward and 
backward filters, the optimum solutions correspond to 
the transfer Function solutions determined in Chapter 2. 
The diflcrence is that the forward filter can now be 
synthesized by one transversal filter with tap spacing 
equal to the Nyquist interval l/B. 

The o p timum so lu t ion de te rmin e d f ro m eqn . (2.37) 
depends on the correlation matrix 


& 

T\. T\ 


“K =K 


'^K S K 

-le-K 


(3.4) 


being non-singular. The quadratic form associated with 
eon. (3.4) is 


QF = u 1 r,, rl u - 2u ’ r,,sA w + w’ s,, s’ 
— K, i\ * k^-k ^ * k • k 


w 


(u» 


T — V ’/ 1 

•-it 




> 0 


(3.5) 


Thus y K is at least non-negativo definite. 

To prove positive definiteness, consider the 
random process 

CO CO 

'K = . S Uii[KT+i/B] - Z w 1 s K _ i 

Irr—oo lor— co 


( 3.6) 
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or equivalently 

*K = a' r K" !l ,S K 

This random process is stationary since both of its 
component processes are stationary. If the power spectrum 
is non-zero in any f requeue/ region for arbitrary | | u| | 2 + 
i | W| | > O, QF is strictly positive and eqn. (3.4) is positive 

definite. The process r(t) has an additive noise component 
with non-zero spectrum over the entire band of interest. 

Also the source correlation in the real problem is finite in 
extent and the source spectrum has energy at all frequencies. 
Thus the QF must be positive unless the process r^ and 
are linearly related. But r^. has a noise component which is 
independent of the source digits, so a linear relation of 
that form cannot exist. Thus Vb y* is positive definite. 

jN. Tin, 

The positive definiteness of V£. not only guarantees 
that the inverse matrix solution to eqn. (2.37) exists, but 
also insures that a gradient algorithm will converge if the 
stop size a is chosen small enough. A gradient algorithm 
for a MMSE receiver takes the form 

f , = jL - a [Wf - \L s' ] (3.7) 

where f is the received vector at the nth iteration and 
-n 

the term in brackets is equal to one half the gradient with 
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respect to f of the MSE expression (2.34). From the 

"nil 

expression for the error process, (3.7) can be cast into 


4 + i = 4-“^ 


n 


(3.8) 


In terms of the forward and backward filter components, 
eqn. (3.8) becomes 


Vl = ^ ~ n ' 3 ' 9) 

fen+1 = 6n + “ 5 e n (3-10) 

The vector r represents the input to the forward filter and 
S represents the input to the backward filter. The compo- 
nents of r and S are the tap values of the respective 
forward and backward transversal filters. Thus the simple 
algorithms (3.9) and (3.10) are implemented by simply 
correlating the tap values with the error process and 
incrementing the corresponding tap gains. From' the analysis 
in Chapter 2 the algorithms will converge if a is chosen 
s o tha t 

J(V V‘ ) < (3.11) 

where J*(.) represents the largest magnitude eigenvalue of 
its argument matrix. In equation (3.8) tap gains arc 
updated after reception of few say L digits. 
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The algorithm (3,8) has been analysed by Monsen [19] 
who determined convergence conditions, the performance 
neighborhood and convergence rate. The performance 
neighbourhood is measured by a quantity M called mis- 
adjustment which is the ratio of the steady state excess 
mean- square error to the IVMSE, i.e. 


M = jT < v n " V 2 > 

where the MMSE is p = l - <VS>* <W l >“ X <VS>. For good 
data performance it is generally necessary to have a small 
misadjustment , i.e, M <<1. 

Mon sen showed that the best choice of averaging 
interval is L* = 1, i.e., no averaging at an. Thus the 
algorithm reduces to 


^n+1 


a \ e n (3.11) 

which interms of forward and backward taps become 


W = a Sn e n 


(3.12) 


— n+1 


A 

b + a S e 
-n % n 


(3.13) 


The forward filter is transversal filter with tap spacing 
equal to l/B where b/ 2 is the channel bandwidth. The 
forward filter output at time nT is given by 
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00 , OO 

Y n = / Y n (t)r(nT+t)dt = ^ Y n [ | ] * [nT+ 


l=-oo 


which can be written as 


(3.14) 


n 


ft ft 


(3.15) 


The forward-filter output has the output of the 
backward filter subtracted from it. The backward filter 
is also a transversal filter with tap spacing T and tap 
gains b^. Its input is the reconstructed data sequence 
which is assumed to approximate the transmitted 
sequence Ki- 

The backward- filter output can be written as 


ft 


y A 

i b i S n-i 


b* s" 
_ 


( 3.16) 


Thus the filter tap gains are updated by simply 
incrementing the present value as a fraction of the product 
of the tap voltage and error sample at that instant. The 
error voltage at time nT is given by 


n 




(3.17) 


The receiver configuration just described is shown 
in Fig. 3.1. In addition to reducing noise and intersymbol 
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Y n = / Y n (t)r(nT+t)dt = |r- z Y n [ \ ] r [nT+ ^ 




(3.14) 


which can be written as 


n 


4 L 


~n 


(3.15) 


The forward-filter output has the output of the 
backward filter subtracted from it. The backward filter 
is also a transversal filter with tap spacing T and tap 
gains b^. Its input is the reconstructed data sequence 
jf which is assumed to approximate the transmitted 


sequence 


The backward-filter output can be written as 


Zn 


X . £ 

i b i S n-i 


/\ 

b ’ 

— ~n 


(3.16) 


Thus the filter tap gains are updated by simply 
incrementing the present value as a fraction of the product 
of the tap voltage and error sample at that instant. The 
error voltage at time nT is given by 


e 


n 




(3.17) 


The receiver Configuration just described is shown 
in Fig. 3.1. In addition to reducing noise and intersymbol 
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interference, the receiver has some additional advantages. 
Since the forward-filter tap spacing is equal to the 
Nyquist interval, the algorithm can adapt to changes in 
transmission delay. Thus, timing jitter introduced by 
changes in transmission delay can be eliminated because 
the delay variations are less than the adaptation rate 
required. 

3.3 DATA RATE: 

The fibre transmission medium has a finite disper- 
sion for tho propagating optical pulse as discussed in 
Chapter 1, and a finite loss rate for the optical trans- 
mitted power. These both limit the data rate that can be 
transmitted as the photo-detector requires a certain mini- 
mum signal power for a given error probability and data 
rate. Semiconductor laser or LED sources have limited 
output power of the order of 2 mw for laser and 100 
for LED that can be coupled into the fibre. These facts 
taken together restricts the fibre length one wishes to 
interface between the source and detector. 

The Fig. 3.2 [21] shows the receiver sensitivity for 
PIN and APD detector digital receivers for various bit 
rates. The figure was calculated for an error rate of the 

o 

order of 10” . The total power available to the system for 




PIN detector under the condition of Fig. 3.2 may be 
described as F(T), and is given approxima tely as 

F(T) =-110 + 4.35t n (T" 1 ) (3.18) 

where 1 T’ is the bit interval. 

This power is shared in various ways in the opera- 
ting system and an operating margin must be maintained for 
the electronics so that small fluctuations do not effect 
the performance of the system. Typically a margin of 5 dB 
is used. The fibre itself is characterized by a loss 
in dB/km, and in addition, there will be power loss due to 
joints, mode coupling, etc. 

Due to the dispersion effect the power transmitted 
in a given bit interval win partly arrive in the corres- 
ponding bit interval at the receiver, and the rest of the 
power will spill over into the nearest neighbour interval 
or oven to more intervals depending upon the impulse response 
of the fibre. This overlap between the adjacent pulses, 
causes what has been referred to as the ISI. Thus the 
overall system can not be designed unless ISI is taken 
into account. The ISI at a finite data rate in the system 
could be reduced by equalization in the receiver after the 
detector, but only at a penalty it pays. The penalty is 
described by F^^ and measured in decibels. Therefore, 
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for a given bit rate, the design equation becomes 
F(T) = 110-4.35 Cr 1 )+P L = a(L) + M+ F ISI 

(3.19) 

where is launched power ( dBm ). VJe will use eqn. (3.19) 
in the next section for link calculation. 

3.4 FIBRE OPTIC DIGITAL LINK CALCULATIONS: 

The Performance of the digital fibre optic link 
depends upon many factors, viz., source, length of fibre, 
attenuation constant of fibre, photodetector and coupling 
loss. As already explained, fibres are dispersive in 
nature, thus pulse propagated through the fibre will get 
spread in time and cause ISI. This will limit the data 
rate that can bo transmitted over a given fibre optic link. 

As discussed in Chapter 1, we have used a PIN detector 
for the system under study. The link calculation for the 
fibre optic communication system under study for different 
lengths and sources are given in Table 3.1. Specifications 
of the components used for link calculation are as follows: 

a) Source; 

i) Laser Diode 

Output power : mw ( 0 dBm) 
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Tabic 3,1: Link Caicula 1; vn for the F.O.C. Digital 
Sys tem 

A. Source-Laser Diode 

i) Fibre length : 2.155 km 

ii) Fibre length : 4.310 km 

Source Power : 1 nW (0 dBm) 

Connector Insertion loss : 2 dB 
Total Power available to system: 

F(T) = -110 + 4.35 -tn (T -1 ) 


Data hate in k bits/ sec 

100 

150 

200 

250 

300 

Total power 

29.87 

28.10 

26.85 

25.88 

25.09 

Connector Insertion loss 

2.0 

2.0 

2.0 

2.0 

2.0 

0 p o ra t in g mn r g in 

5.0 

5.0 

5.0 

5.0 

5.0 

Power available to 
trnn smis s ion me dia 

22.87 

21.10 

19.85 

18.38 

13.09 

Fibre loss for 2.155 km 
(4.5 dB/km) 

9.69 

9.69 

9.69 

9.69 

9.69 

Power available for 

I SI penalty ( dB) 

Fibre loss for 4.310 km 
(4.5 dB/km) 

13.18 

11.41 

10.16 

9.19 

8.40 

19.39 

19.39 

19.39 

19.39 

19.39 

Power available for ISI 
penalty (dB) 

3.48 

1.71 

0.46 

- 

- 


B. Source-L^D 

i) Fibre: .length : 2.155 km 

Source Power : 50 p'«/( -13.01 dBm) 
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Total power available to 
Date Ratos in M bit^ 
Total power 

Coupled power to system 
Connector insertion los 
Operating margin 

Power available to 
transmission media 

Fibre loss for 2.15b km 

Power available for ISI 
penalty ( dl3) 


110+4. 35^ T -1 ) 


75 

100 

125 

31.12 

29.87 

28.89 

13.01 

13.01 

13.01 

2 b 

2 o 

2 *'•> 

5 • o 

5 -o 

5-o 

ii.ii 

9.86 

8.88 

9.69 

9.69 

9.69 

1.42 

0.17 



system s F(T)= - 
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ii) LED - IRE - 161 ( LASER- DIODE) 

Output Power : 50 pw at A A = .2 

Current : 100 mA 

Core : 50 pm 

b) Fibre - (M3 -05): 


Length 
Attenua tion 
NA 


2.155 km 
4.5 dB/km 
0.23 


Core 


63 pm 


Material dispersion 


90 p.sec/km-nm 


c) Connector 0FP-101 


Insertion loss = 1 dB 

d) Detector 

PIN Diode - 0FQP-04 (ITT, UK) 

Dark current : 0.2 nA 

Sensitivity : As shown in Fig. 3.2 for 

BER 10“ 8 

Responsivity : 0.5 A/w 


It is assumed that the connector insertion loss is 
same for the laser diode and the LED. Also connector loss 
for the two sections of fibre is assumed to be negligible. 
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From the link calculations we observe that with the 

laser diode the system can work up to a data rate of 

200 M-bit, for a fibre length of 4.310 km whereas with the 

ra te 

LED the system can work only upto a data/of 100 M-bits 
for a fibre length of 2.155 km. As for data rates higher 
than this there is no power available for meeting the ISI 
penalty. Therefore, to achieve tho given BER of 10 -8 at 
higher data rates, either one should use the laser diode or 
decrease the length of tho fibre. Link calculations with a 
laser source for a fibre length of 2.155 km varifies the 
latter suggestion. 

3.5 IMPULSE RESPONSE: 

Impulse responses of the various parts of the fibre 
optical system and the overall link have been given in the 
following sections. 

3.5.1 General: 

Impulse response of an optical fibre is generally 
found by measurements of the delay distortion. In the 
system under discussion a graded index fibre is used, 
and the measured impulse response of the fibre provided by 
the manufacturer is shown in Fig. 3.3. It is assumed that 
the impulse response is normalized and the time is given 
in nano-secs. The suppliers have not given the exact value 

Mo, f 
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of the magnitude (scaling factor). The value of the 
scaling factor is found as follows? 
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We know - 

oo 

/ (Kh) r dt = 1 A* small time interval 
o 

or, 

K 2 Z h 2 i = 1 (3.20) 

For calculating K,.\_ is taken asO.l ns. and the value 

q 

from calculation turn ry'- to be 1.5b x 10 . As we have 
taken only a finite number of samples of the impulse response 

q 

the exact value of K is 1.5x10 . 

3.5.2 Impulse Response of the Electrical System: 

The optical signals are converted to electrical 
signals by the PIN diode detector and then these are ampli- 
fied by the pre-amplifier, as the data rates of interest 
are high, a high speed FET ( CA 3127E) is used for the pre- 
amplifier circuit. Tho characturis tic of CA312.7 and the 
PIN diode are as under - 

CA3127E 

Input resistance = 400 Si. 

Input capacitance = 3.7 pF 

PIN Diodo 


Ca pacitance 


1 pF 
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The equivalent electrical circuit of the PIN diode and the 
pre-amplifier is given in Fig. 3.4. 



£ i q ■ 6. 4- 


Current due to the shot noise and dark current are 
not considered being negligible. The bias resistance of the 
photodetector is high and shunted by the 400 SL. input 
resistance of the pre-amplifier. The resultant values of R 
and C of the model are 

R = 400 -0. 

C = 4.7 pf 
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fhe impulse response of the detector and the pre- 
amplifier is then 


H(s) 


V 2 (s) 

vJiT 



1 

c 

S+ 1/R C 


(3.21) 


;.h(t) = 1/C . e” l/R c t (3.22) 

Impulse response of the equivalent circuit is plotted as 
shown in Fig. 3.5. 


3.5.3 Impulse Response of the System: 

To evaluate the performance of the DFE, the total 
Impulse response of the system is considered. The fibre 
optical communication linxs under consideration are as 
follows : 


a) Source-Laser Diode, Detector -PIN Diode Link length = 2.155 km 

b) Sourcc-LED, Do tec tor~P IN Diode, link length = 2.155 km 

c) Source-Laser Diode, Do tec tor -PIN Diode Link length = 4.310 km 


a ) Laser-piN Combinati on. Length =■ _2.,1.55_km: 

The material dispersion of the graded index fibres 
is 90 p. secs/km. -nm. It is assumed that the spectral 
width of the laser is 2 nm. Thus the total material dis- 
persion duo to the laser is very small, i,o,,0.45 nsecs. 
Therefore, the effect of material dispersion due to the 
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source in this case is neglected while calculating the 
total impulse response of the system. Vie know the impulse 
response of the fibre and the electrical circuit (detector + 
pre-amplifier). The total impulse response is found by 
convolving the two impulse responses. The total normalized 
impulse response is shown in Fig. 3.6. 

b ) LED- P1N Co mbinatio n Length = 2.JL56 km : 

Optical fibre properties depend upon the wavelength^ 
of a light source, i.e., g(t) = g(t,X)> where g(t) is 
the impulse response of the fibre. The spectral width of 
the LED is 40 nm. Therefore, the material dispersion will 
be much more compared to the laser source. The LED, IRE-161 
of Laser Diode Lab. is used in the fibre optic communication 
link under discussion. The output spectrum of the light 
source is as shown in Fig. 3.7. 

The resultant pulse response of the fibre link 

is [22] 

g r (t) = / p(X)g(t, A)dA (3.23) 

* *-c© 

whore p(A ) is high source spectrum. 

The resultant impulse response is shown in Fig. 3.8. If 
p(A) is not available for a LED, wo can use equation (3.24) 
whore a Gaussian intensity spectrum for the LED has been 

assumed. 




Result ant impulse £ esPon sc~. 
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WAVE i_£ V <$Tt4 w m 


F/c> 3-7 


OUTPUT ZPBCTRA 




resultant impulse Response. 



G m (f) = exp I - | . (2itfML a) 2 
or y 

^ “ f2ixMLa*“ exp j^“ \ (t/MLc) 2 J (3.24) 

v^he re M = material dispersion in ps/km-nm 

L = length of fibre in km 
a a spectral width of LED in nm 

The overall response win then be obtained by convolving 
9 m ( t) with g( t) . The output spectra of the light source 
IRE - 161 is not Gaussian. The value of a is found by 
curve fitting. It is 14 n.secs on the right hand side of 
the peak value and 17 n.sec. on the left hand side. In 
our case the resultant impulse response of the link is 
obtained numerically by using the actual p(A) as provided 
by the manufacturer. 

The plots of Fig. 3.9 and 3.10 have beefr obtained 
using the Gaussian assumption and using values t)f a of 14 and 
17, respectively. The resultant impulse response shown in 
Fig, 3.8 corresponds to a worst case situation. 

c ) Laser~PIN Combination. Length- 4.310 km : 

The total impulse response of the system is found by 
convolving the impulse response of the fibre of 2.155 km 
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lenqth with itself and then with the electrical circuit 
impulse response. Tote.] : i ' 1 Ls e response is shown in 
Ficj. 3 » 1 ,1. which is used co find the Probability of error 
of the system of fibre lenqth 4.310 km. 

3.6 CHANNEL BAND-WDTH OF THE SYSTEM: 

As discussed in Chapter 1, the system can be divided 
into two parts, i.o, optical and electrical. 

The* system under discussion have electrical as well 
an optical portions. The 3 dB channel bandwidth of the 
combined system is calculated. 

It is observed from Section 3.2 if the channel 
is band limited to B/2, it is reasonable to position a filter 
prior to the equalizer which will also band limit the 
addi U vu noise to the range |fj < B/2. Therefore one should 
know the channel bundwid'.; of the system upto the input to 
the equalizer. 

Vic have* the pre-amplifier before the equalizer, so the 
additive noise is also limited to tho bandwidth of the pre- 
amplifier. Therefore, an additional filter to band limit 
the additive noise is not required. 

The impulse response of the fibre or the overall 
system in our case is known. For the impulse response 
h(t), 0 i t t the transfer function is 




C u h n.<vh o S) <t c . 

‘Resultant /mpulsb response 
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r y 

H(f) = / h(t) cos 27rf tdt+- j / h( t)s in2itf tdt 

0 o 

= A(f) + j$(f) 

| H ( f ) J 2 = A(f) 2 + B(f) 2 

where A(f) = f h( t) cos2itf tdt and 
o 

y* 

B(f) = / h( t)s in27tf tdt 

o 

The integral functions are computed numerically. The 
3 dB channel bandwidth is obtained from the plot of 
i H( f ) | 2 * 

The magnitude of the transfer functions of the 
electrical circuit, the fibre optics, and the total system 
(a) are plotted in Fig. 3.12 and Fig, 3.13. The values of 
the 3 dB bandwidths are obtained as under: 

i) Electrical Circuit = 85 MHz 
Fibre optic = 45 MHz 

Total bandwidth of system (a) = 38 MHz 
Total bandwidth of system (b) =24 MHz 

3.7 DESIGN OF DFE: 

In the following the choice of the step size (a) and 

the tap gain coefficients of the DFE have been dealt with 

% 

in brief. 




kn JjJi nrn OF Pt vF-amdi tptco JZ 
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3.7.1 Effect of Step Size (a): 

The critical consideration in the choice of the 
algorithm parameters is to maintain a misadjus tment small 
relative to unity. 

The algorithm parameters which can be varied are 
the step size a and the averaging interval. We have 
already assumed the averaging interval to be unity, a 
should make the receiver tap gains move to the optimum 
neighbourhood from an initial setting corresponding to the 
absence of an equalizer. Thus the problem is to maximise 
the convergence rate. 

While carrying out simulation on the digital computer, 
it is observed that if we make a larger then the probability 
of error is more, though convergence is faster, on the other 
hand if the step size is kept too small then the convergence 
rate is very slow. The algorithm will converge if a is 
chosen so that 

< | 

D.A. George [23] has used a value of a as low as 0.016. 

In this thesis we assumed that the reconstructed data 
sequence has no error. Therefore, a should bo so chosen 
that the probability of error is low. When the equalizers 
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tap gains can adaptively move to the optimum neighbourhood 
from the initial setting, the equalizer is said to have 
achieved lock. 

In the simulation the tap gains of the forward and 
the backward filters are initially set to zero except for 
the tap corresponding to F Q (main tap in the forward 
transversal filter), which is set to unity. It is noticed 
during the simulation that if after a number of iterations 
the weight of F Q is around unity, then the best tap gain 
settings are achieved. But in case the value of F Q tap 
converges to zero then lock is not achieved. The error rate 
becomes very high after a few iterations. The forward 
filter tap weights sot to zero and backward filter tap 
weights set to such a value that their sum tends to bo 1 
or -1. In this case the reconstructed data is independent 
but the error, o^ = ^ is identically zero for all 

times. Therefore, the selection of a is very important and 
should be chosen so as to achieve the lock condition. Also 
the tap weights of F q should not converge to zero. The 
value of F q should bo around unity to get the best results. 

3.7.2 Tap Gain Setting: 

Length of transversal filter : The parameters of in teres 
in the DFE are the number of taps in the forward filter, i.o. 
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~N to N, the forward filter tap spacing v-p and the number of 
taps in the backward filter, say M. The optimum DFE 
requires N = oo = M and Y = ^ , where l/B is the Nyquist 
interval. In practice the transversal filter will be of 
finite length. 

Monsen [24] in his studies kept tap spacing for the 
forward filter as 66 ns (corresponding to 15 MHz), data 
rate between 1.5 to 12.6 M bits/ sec. and three taps in 
the forward filters. In the system under study the channel 
bandwidth is about 38 MHz, and data rates much higher than 
the channel bandwidths are used in simulation. From the 
simulation studios, the following observations can be made 
about the- number of taps and their spacing. 

i) When the forward filter taps are spaced at the symbol 
interval, the probability of error is loss than or 
similar to the value obtained when the taps are spaced 
at the Nyquist interval. 

ii) If we keep the number of the forward taps equal to the 
number of future interfering samples plus one then the 
probability of error is smaller. 

iii) Any additional tap added in the forward filter after 
the main tap F 0 does not reduce the probability of 
. The results are tabulated in Table 3.2. 
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CHAPTER 4 


SIMULATION OF DFE AND RESULTS 


4.1 GENERAL: 

To evaluate the performance of DFE, the actual system 
is simulated on the DEC-1090 system. The overall impulse 
responses of the link are considered for the simulation of 
the system. The simulation consists of finding the tap 
gain weights for the forward and backward transversal filters 
and then using these tap gain weights to find the probability 
of error at different data rates for various SNR values. The 
probability of error is determined for 2.155 km and 4.310 km 
link lengths for SNR values from 6 dB to 14 dB. Performance 
of the DFE for these link lengths is evaluated considering 
the different optical sources, i.e., laser — diode and LED. 
Maximum data rates possible for a certain probability of error 
without an equalizer have also been found. 

For simulation random numbers with uniform distri- 
bution (data) and Gaussian distribution with zero mean 
(noise) arc generated by NAG sub-routine. The uniformly 
distributed random numbers are convolved with the total 
impulse response of the system, and the Gaussian noise after 
convolving with the pro-amplifier and detector impulse 
response is added to it. Now the signal is fed to the DFE. 
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Signal flow diagram for the simulation is shown in Fig. 4.1. 
Error signal is used to update the tap weights during the 
adaptation process. Once the tap weights of the DFE are 
fixed, error signal e is not used* 

4.2 SIMULATION RESULTS: 

4.2.1 Source-Laser, Fiber Length = 2.155 km: 

The simulation is carried out in two parts. First 
we find the optimum tap weights for the forward and backward 
filters for the different data rates at various SNR values 
by the adaptive method. As discussed in Chapter 3, the 
forward filter having number of taps equal to the future 
interference samples plus one and spaced at the symbol 
interval is considered in evaluating the optimum tap weights 
and the probability of error. For various values of SNR 
the noise variance and the standard deviation is given 
in Table 4.1. 

The optimum tap weights are determined by carrying 
out a number of simulations with different values of step 
size and then freezing the values of the tap weights at 
different iteration numbers. The tap weights which corres- 
pond to the minimum number of errors in the above simulation 
is used to find the probability of error of the system. The 
results from one such simulation are given in Table 4.2, 
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Table 4,1; No Ice Variance and Standard Deviation 
at Different SNA Values 


SNR 


Noise Variance 


Standard Deviation 


6 

0.2512 

0.5012 

8 

0.1585 

0.3981 

10 

0.1 

0.316227 

12 

0.0 631 

0 . 251188 

14 

0.0398 

0.199526 

16 

0.02512 

0.15849 

2d 

0.00398 

0.0 631 

32 

0 .000 63 

0.02512 
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Table 4.2 

SNR =10 No. of sample processed = 5000 
Data Rate 150 M bit/sec. 


Step size Error Freezed Tap weight ______________ 

value at Forward Backward 

Lli O) LZ 1 (-21 


0.10 

7 

495 

0.9803 

0.08568 

-0.0 6226 0.07677 

0.1 

8 

990 

0.9741 

0.004077-0.06914 0.1162 

0.1 

9 

1500 

0 .8070 

-0.01352 

0.10 62 0.1972 

0.05 

4 

495 

0.9510 

-0.04138 

-0.01124 0.1553 

0.05 

5 

65 

0.90 28 

-0.003379 0.02807 0.18607 

0.05 

6 

95 

0.9244 

0.05050 

0.007390 0.1951 

0.65 

6 

145 

0.8819 

0.03280 

0.06871 0.1385 

0.05 

10 

195 

0.9123 

-0.06046 

0.06989 0.1145 

0.05 

8 

295 

0 .8847 

0.04733 

0.03952 0.1671 

0.1 

9 

195 

0.9410 

-0.0 6926 

0.04480 0.1248 
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The optimum tap weights for the forward and the 
backward filter for the different data rates at various 
SNR values are given in Table 4.3. 

Now the DFE is synthesized with the tap weights as 
calculated above. The average probability of error at data 
rates from 150 Mbps to 375 Mbps and for SNR values from 6 dB 
to 14 dB is determined. These are shown in Table 4.4 and 
plotted in Fig, 4.2* It is evident from the curve that 

n 

for 2*155 km fibre a probability of error less than 10 can 
be achieved at higher SNR values. Table 4.4 also shows the 
total number of symbols processed. As discussed in Chapter 2 
the error propagation effect is less as the ;ENR increases. 

In Fig. 4.3 we have reproduced the error occurance during 
the simulation for SNR values of 6,8,10 and 12 dB. These 
are a sequence of ' and 1 ** . 1 .* indicates a correct 

decision and * indicates an error. We notice that error 
bursts are more at higher data rates and lower values of 
£^R. If a group of errors is followed by correct decisions 
and the number of errors is equal to the number of backward 
filter taps then a particular error burst can be considered 
to be over. Any subsequent errors are duo to noise and I SI 
arising from the future samples. 

4.2.2 Source-LED, Fibre Length = 2.155 km: 

In this case the source used is LED. • Therefore the 



Table 4.3 : Tap weights for source-Laser, link length =2.155 km 
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Table 4.4: System Performance (Source-Laser, link 
length = 2.155 km) 


Data 

Rate 

M bit 

SNR 

(dB) 

No. of symbol 
processed 

Error 
due to 
ISI only 

Error 
wi tho ut 
equalizer 

Error 

P(e) 

150 

6 

50 ,000 

0 

1542 

1202 

2.40 xl0“ 2 


8 

50,000 

0 

551 

313 

6. 26x10” 3 


10 

10 5 

0 

269 

86 

8. 6xlO” 4 


12 

2x10 5 

0 

72 

3 

1.5x10” 5 


14 

8x10 6 

0 

168 

0 

Better thai 
1.25xl0" 7 

225 

6 

50,000 

0 

2812 

1271 

2. 54xl0” 2 


8 

50,000 

0 

1645 

390 

7 . 8xlO” 3 


10 

10 5 

0 

1532 

10 6 

1.0 6x10" 3 


12 

10 5 

0 

564 

3 

3x10” 5 


14 

10 6 

0 

1542 

1 

1x10" 6 

300 

6 

50,000 

r, 

4836 

1954 

3.9xl0" 2 


8 

50,000 

0 

3627 

519 

1.0 38x10” 2 


10 

5 

10 

0 

5304 

130 

1.3x10” 3 


12 

10 5 

0 

3780 

13 

1.3xl0 -4 


14 

io 6 

0 

17511 

11 

l.lxlO” 5 

375 

12 

10 5 

5536 

8792 

21 

2. lxio” 4 


14 

10 5 

5460 

7977 

2 

2 x 10" 5 
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material dispersion due to the source is much more and it 
spreads the impulse response further in time domain. The 
resultant impulse response of the system is shown in Fig, 

In this case also I i.rat we determine the optimum 
tap weights for SNR values of 1Q,12 and 14 dB, The optimum 
tap weights are given in Table 4.5, Now the probability of 
error is found by using the tap gain values as determined 
above. Probability of error for a^R values 10,12 and 14 dB 
are plotted in Fig, 4,4, It is clear from the figure that 
the probability of error loss than 10“ 7 can be achieved at 
higher SNRs and lower data rates. Table 4,6 shows the total 
number of symbols processed, the number of symbols received 
in error and the probability of error. 

In Fig, 4.5 we have reproduced the system output 
photograph with errors. These are again a sequence of 
and , * t . We observe that the error propagation effect due 
to backward filter is very loss. 

For the system under discussion we have also carried 
out simulation at reduced data rates for 3\IR values of 10,12 
and 14 dB, when the equalizer has not been used. In this 
case, the output of the pre-amplifier is fed directly into 
the decision circuit. It is seen that at very low data rates 
wo can dispense with the equalizer as the symbol interval is 
very large and there is much loss ISI.t Table 4.7 shows the 



Table 4*5: Tap weights for Source-LED* link length « 2*155 km 
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Table 4.6: System Performance (Souxce-LED, link 
length = 2.155 km) 
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Table 4.7: System Performance without DFE (Source-LED, 
link length = 2.155 km) 


Data Rate 
in M Bits/s 

SNR 

(dB) 

No. of Symbol 
pro cessed 

mm *m. i *r m vmrwm 

Error 

P(e) 

75 

10 

10,000 

32 

3. 2x10“ 3 


12 

50,000 

9 

1.8x10“ 4 


14 

5,10, 7; 5 

5 

9. 8x10” 6 

62.5 

10 

10,000 

9 

9xl0“ 4 


12 

5 

10 

10 

IxlO” 4 


14 

4x10 6 

5 

1.25x10" 6 

50 

10 

20,000 

15 

7.5xl0“ 4 


12 

111,762 

5 

4.47xl0" 5 


14 

5x10 6 

1 

2xl0" 7 



93 


data rate, number of symbols processed and the probability 
of error without the equalizer. 

4.2*3 Source-Laser, Fibre Length = 4.310 km: 

In the system under discussion we assumed two equal 
fibre sections of length 2.155 km. each have been ideally 
connected to one another. For simulation the impulse res- 
ponse of the system as shown in Fig. 3.11 is considered. 

Again for this system first we determined the optimum tap 
weights of the forward and backward filters by the adaptive 
method for SNR values of 10,12 and 14 dB. The optimum tap 
weights are given in Table 4,8. 

The average probability of error at different data 
rates for 4.310 km fibre length and for various SNR values is 
determined. The probability of error for various 3MRs is 
plotted in Fig. 4.6. It is observed that with the increase 
of the link length the probability of error for the same data 
rate and SNR is higher. Table 4.9 shows the number of symbols 
processed, the number of errors in simulation and tho proba- 
bility of error. 

For this system also we have carried out simulation at 
much reduced data rates when tho equalizer is'*not used. Tho 
results are shown in Table 4.10. The main interference is 
due to tho noise. Therefore, tho system works bettor without 
any equalizer at higher values of SNR. 


Table 4,8: Tap weights for Source-Laser, link length = 4,310 km 
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Table 4.9: System Performance (Source-Laser, link 
length = 4.310 km) 


Data Rate SNR No. of Error Error Error p(e) 

in (dB) symbol due to without 

M Bits/s processed ISI equalizer 


only 


150 

10 

5 

10' 

0 


12 

2x10 5 

0 


14 

5x10 6 

0 

200 

10 

10 5 

0 


12 

10 5 

0 


14 

2x10 6 

0 

250 

10 

50,000 

4690 


12 

50,000 

4585 


14 

10 6 

93772 


1704 

109 

1.0 9x10“ 3 

1294 

5 

2. 5x10“ 5 

9355 

3 

6x10” 7 

6428 

121 

1.21x10“ 3 

4934 

6 

6x10“ 5 

76110 

3 

1.5x10“ 6 

5640 

152 

3. 04x10“ 3 

5240 

7 

1.4x10“ 4 

102037 

11 

• 1.1x10“ 5 
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Table 4.10: System Performance Without DFE 

( Source-Laser, link length = 4.310 km) 


Da ta Ra te 
in M bits/ s 

SNR 

(dB) 

No. of Symbol 
Pro cessed 

Err 

P(e) 

100 

10 

10,000 


6.9xl0" 3 


12 

10,000 

14 

1.4x10"* 3 


14 

10 5 

9 

9x10" 5 

87 

10 

20 00Q 

43 

2. 15x10" 3 


12 

2.0 000 

9 

4. 5xl0" 4 


14 

10 5 

8 

8x10" 5 

75 


20,000 

31 

1.55xlO" 3 


.2 

10 5 

16 

1.6xlO" 4 


14 

10 6 

2 

2x10" 6 

60 

10 

20,000 

■30 

1.5x10" 3 


12 

46,300 

5 

1.08X10" 4 


14 

10 6 

0 

Better than 10 


6 



CHAPTER 5 
CONCLUSION 


Ibis thesis has been directed to evaluate the per- 
ioriihnce of some fibre optic digital communication systems. 
io have shown an optical fibre system with its essential 
eomponen ts. The choice of source and pho 'to detector depends 
upon the length of the fibre, attenuation of the fibre and 
desired probability of error for a given data rate. We 
have soon that to uso the same fibre length at higher data 
rates, one has to use a laser diode as the source, as the 
power required to combat the ISI penalty is not sufficient 
when we uso a LED, Similarly the increase in fibre length 
limits the data rate which can bo transmitted over the link 
for a desired P(e). 

We have considered systems with different sources 
and fibre lengths as follows. 

a) Source-Laser diode, link length = 2.155 km 

b) Sour ce-LED, link length = 2.155 km 

c) Source-Laser diode, link length = 4.310 km 

\k. have noticed that the fibre channel spreads the 
optical signal in time domain. This dispersion and additive 
noise of the system cause Id* The effect of ISI can bo 
reduced with an equalizer. 



99 


In this thesis we have designed a DFE based on the 
Mon sen * s Gradient-Technique. DFE uses both the received 
signal and reconstructed data to form its decision. The 
tap weight values for the DFE have been determined by 
Monsen's adoptive linear filtering method. 

To determine the tap weight values and to evaluate 
the performance of the system, we carried out simulation 
on the digital conputer DEC-1090 system. Optimum tap 
values are obtained when the number of taps of the forward 
filter is equal to the number of future ISI sample plus 
one and taps are spaced at the symbol interval. 

The initial aim was to find the maximum data rate 
and SNR value for a system which gives a p(e) of 10“ 8 . To 
get the P(e) of 10~ 8 , one has to generate and process on 

Q 

an average 10 noise and data symbols. From simulation 

fs 

wo have noticed that to generate and process 10 symbols 
on the conputer CPU and execution time are approximately 
6,8 and 35 minutes, respectively. Therefore, a maximum 
of 8x10^ symbols could be generated and processed in the 
simulation. 

The fundamental assumption in the design of the 
UFE is that the decisions arc error free, but, in every 
simulation certain amount of error, i.e. 'Error Bursts* 
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have o ccurcd. Any wrong decision affects the decision 
i.n the DFE till it propoi 'os in the feedback filter. 

Error hursts are more at higher data rates and lower SNR 

values. 

From the results we have noticed that for a data rate 

of ISO in bits/ sec., the system with the LED source gives 

a probability of error of 2.5xl0~ 5 for a SNR of 14 dB. 

Whorc-a s, with the laser diode a probability of error of 
-7 

1,25x10 is obtained for the same SNR value. Also from 
>-hu extension of the performance curve (Fig. 4.2) for 
150 M bits/scc. it is observed that a P(e)of 10 for the 
fibre length of 2.155 km can bo achieved at an approximate 
SNR of 16 dB. Whereas for the fibre length of 4. 3lO km 
(Fig. 4,6), it can bo achieved at approximately 18 dB. 

For systems (b) and (c) performance of the links have 
also been evaluated without an equalizer. For source~LED 
and a link length of 2.155 km, it is observed that to 
achieve the same P(e)at a SNR of 14 dB, the maximum data 
rate which can bo transmitted without an equalizer is 
npproxinva tely &>A of the data rate obtained with the 
equalizer (Ref. Table 4.6 and 4.7). 

For the source-laser and link length of 4.31 km, 
to achieve the same P(e) at a 34 R of 14 dB, the maximum 
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cint*'. im to which can be transmitted without an equalizer 
* Igog than 50^ of the data rate obtained with the 
equalizer (Ref. Table 4.9 and 4.10). From this we 
o b.nerve that with the equalizer, we can transmit atleast 
double the data rates than without the equalizer to 
achieve the same P(e). 

Capt. R, Wahi [15] designed a DFE based on Austin’s 
sub-optimal method. He .wa used a data rate of 10 M-bits/sec. 
and a fibre whose inpulse response at 10?< amplitude has 
a full width of 275 ns. For 14 dB SNR he has obtained a 
P(e) of 3xlO~~ > . In our study we have used a resultant 
inpulse response, whose width at 10 /• amplitude is 10.5 ns. 

We have used data rates from 150 to 375 M-bits 

In our case for a SNR of 14 dB, the P(e) 

at 300 and 375 M-bits/sec. is l,5xl0~ 5 and 2xl0“ 5 , respe- 
ctively. For 150 M- bit/sec. and SNR of 14 dB, we get P(e) 
better than 1.25xl0~ 7 . In Capt. R. Wahi’ s study the ratio 
of full width at 1 Oyi to symbol interval is 2.75 and in 
our study it is 3.18 and 3.9 for data rates 300 and 
375 i,l bits/sec, respectively. From the above this system 
scums to be better. 

One could study the performance of DFE after actual 
measurements of the overall inpulse response of the system 
and then simulating on the computer. 
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Another issue for further study is the simulation 
of the entire system including line coder, bit synchronizer 
and decoder etc., on the digital computer. The performance 
of the system can also be studied after hardware imple- 
mentation of the DFE and then its performance can be 
compared with the simulation results. 

Finally system can also be investigated by considering 
tho affect of shot noise on the system. 
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START 


Read initial data, (No. of forward and 
feedback taps, set F Q unity, step size 

SNR, counter ' 


Generate Data 


Convolve data symbol with samples 
of system impulse response 


Generate noise sample 


Sot iteration No. where 
adaptation stop 


Compute signal through 
adaptive DFE 


Compare with Input 


Increment error 
count 
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t’l', "ofU'i A ( 1 0 «: ) , FF (SO) »R|jf‘) 0 ) ,[')( tOOOuO) 

in a 1 1 a.ik 


■ , ” 1 ‘ 'i* T »' = IS * , K T IjFc * T A P 2 . 0 A T ' ) 

r n«i"' ,( ; 'rTst9,l>tftfTCK=*DSN # ,FU.K='T.\P22.0AT') 

= K^kO-I GO Mti7 

nnj ;■'*( i*j ,<•>:») 

Fit'' , A I T 2U a. , # .9f)Rs 1 0 FRli^st SO Vh? AT. PA (STEP STZF) =U .OS # ,/20X 

Af.f* * = ,s'f ( P TtZf, 

M, 1 ’ \S', ,t>S 
ft* vu'i t.ti 

^TPUAs.'OfA JOAPti DEVIATION (FOR flNR) 

r,i«* 'A=".toPi>?b 

P( T japPh-AUpr.lFIFR IMPULSE RFSPONSF SAMPLE V ALOES 
r *’( 1 1 = 1 . 0 00 ? P ( 2 ) =0 .0339 

^ms'1.0 ? Go (2) =0.0 J«N(3) = 0.0 ? GN (4 )*0 .0 7 GW ( 5 )*Q . 0 

’UTjsjYtfTfcM I MPUIjK RFSPUNSF SAMPLE VALUES 

>»( 13=0.09 H? 0(2 ) = l, Of H(3) = 0.S007J HC4)S0. 3123, ul(5)=0. 2102 

u(*)=O.OO27;U(7)s0,0190 

Xt U«»* .ur Xf 2)s0,O;X(3)s0.0lX(4)*0.O* JT(5)*0,0;X(6)s0.OjXC7)s0.0 
*(S)a',...>|XmaO.OjXUO)»0.njX(U)eO,0»X(t2)*0.0;Tf(l3)sO.O 
T ( T ) sijHC A f TU V OF SIGNAL IN FORWARD FILTER 

T( U=o.tt? r(2)aO.OOOO? T(3)=0.0000;T(4) = 0.0000;T(5)=0.0000 

T(S J8O.on(M)#n.0j T( 10)80.0 
F(T)*pnaWARD FILTER TAP VALUE 

F ( 1 3=1 .0000? F( 23= .0000 0?FC 33=0.0000000? F(4)=0.0?F(5)*0.0 

Fl7)*0.0?F(8)=0.0?Ft9)*0.0:F(10jry.O 

S(n*LnCATIOM OF BACKWARD FILTER FOR RECONSTRUCTED DATA 

S(t 3*0.0? SC2)»0.0 ?$C3)*0.0 ?S(4)=0.0;S(5)*0.0?S(6)=0.0?SC7)*O.Q 

S(ft)*0.0?S(9)*0*0?S(10)=0.0?S(ll)=0.0 

Rf I)«FEEDBACK FILTER TAP VALUE 

R (1 3*0 .OO0QO0J 0(23=0.0000 ?B (3) «0. 000000? R(4)*0.0000fBC5)»0. 0000 

R(6)*0,0 10(73*0.0000 
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,7'>= ;r> OK SAMPLES OF PRE-AMPLTFTK9 IMPULSE RESPONSE 

Jii = ‘ " i'M'JP OF oAMPl^G l)F PYSrFm IMPULSE RESPONSE 

-IF = FIin|.;R OK FORWARD FILTER TAPS 
Ti>= jIMP.KR H K FFKimAC* FITTER TBPS 
ThSsNiF’.trK UP FUTURE ISi SAMPLE 

T FTP w af T 0 w ohr up future isi SAMPLE REPORE MAIN tap of FWD r tL v EP 

r-r 1 T=iP,.,+ £rxR'., 

IHA 10i- -ITER.Vrin-j UO FROM fl/FEPE VALUES OF TAP GAIN STORED 

TLA RVslTF'MTin i Jo URTU *”irn VALUES OF TAF GAIN SIGHED 

fFT.ia ITERATION Jo aT ‘nHTEH ADAPT A TION Sl'Op 

.u»=,i? 113/ ; UF=2; jn=S; IES=1 ; if rOM = l ?TTPT=2 

TUA.iC'-s 14 i n ; URAMGEa 1497 ; IFIX= 149S 

Mr 1 Oil'jit 

Nl'SVf 

VsL 

UK 1=0 

MOj.San.0 

Mn?sU .0 
r Ali 1 1 GUSCCF 

no to I = l,'J 

KrFti 

AKSGUOF AFU) 

TFCAK.Oe. 0.5) GO TO 31 
AK=-t ,<» 
c;n i-n vi 
AK = t .»> 

nil ?i >T=1 , ITOT 

A(.T) = riFUT DATA TO THE SYSTEM 

A(.T)sAf J + l) 

rot? TIM HE 

ACT COf +1 ) = AK 

PUT JsUUMRER ONE LESS THAN THE 610 OF H SAMPLE 
00 33 

X(J)«X(J + n 

CONTINUE 

XCJtO*AK 
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h - 1 ' » 1 1 

r - "Hi ilK ’ r i’P'lTin* iCOij A |j TO HU OF u SAMPLE 
1)0 34 F = 1 ,,Ui 

TKCt r- u.or.o J on rn 20 

VKs;Kfi;i U *x(,m + W) 

3 4 r* i r ? 1 ?^ 

r F JO A to K ri OF FUTURE 3 A 1PLE III COMBINED TMP RESPONSE 

°9 TKfK.M’.IPS) GO TO ID 

GK^GubnuFC MiE/vH,SIG,JA) 

<’ .T r.i it Hi‘0 Mf) t IF PRO A3P SAMPLE MINUS ONE 

l>0 3’> 1 = 1 , JP-1 

r-iCJ-UF ,(Jf 1) 

35 r.pi fTM'i ; 

G,C.JP)= ;k 
AhF=h . n 

r mu of rjy.tvrruF k-juaTi to no of pre-amp sample 

DO 3b T=1 , TP 

TF(CT-TFS-.T).LT.O) GO TO 30 
ANK = A tlK tR CtT ) *GH ( Jp + 1 - J ) 

3b CUHTTO’IE 

30 TK( ^F.GIJ.O) GO TO 60 

FK=-1 .0 
GO TO 61 

60 FK = 1.«' 

61 TFU-f 1 + IFTIUD.E0.EK) GO TO 62 
,TE1 s.Tol +1 

C PKsINPHT otGMAL TU OFF 

62 RK s^K+ANK 
TEfBK.GL.O) GO TO 3? 

C SK = !<ECO JSTRUCTED DATA WITHOUT DFE 

SK =-1 .0 
GO TO 38 
37 SK =1 .0 

3b TF(A(1+IFTRM>.EQ.SK ) GO TO 40 

N()13»N01S+t 

TF(K.GT.tO) GO TO 41 


40 



JP T j'K t "i >j , 9‘! } C VK, AWK,UK,Xl JH> ,GN (JP ) , A ( l +IFTRM) , SK , AK , OK ) 
‘‘0 rf t'K, 3n;t t j.5),*>X,2FH.3,2F5..1,5X,Ai,F5.1 .2X.F6.3) 

r tap wu iteration start kpom this point 


•1 3 on t 7 r=l f ,117-1 

TC T)=TC,J+n 

12 f ii H TT 

Tf.TKisoi*; 

C KuUATl- ,s TO uu OF FUTURE SAMPLES +NU UF FWD TAPS BEFORE **ATN TAP 

TFfrv.LF. f tTUTJ ) 00 TO 10 

r orvt. HAoc.i: ri which vat.vf of tap gain reoutrfd 

TFO .uF.TOA.'iOK) 00 TO ?h 
TFIK.i,”'. TRAdGC) GO TO ?h 
no 11 ,T=1,,1F 

rpfJ) =FC-T) 

11 roOTlfinc 

f.»U 0 0 

R0(j1 SB(,7) 

OJ C()NTTm"K 

r YRsrulfTPtJT OF FORWARD FILTER 

?fo YK=o 

00 1 3 fat , IF 

Yh*YKtTU)*Ff J) 

13 rijNTTfini-: 

on 14 T=S 1 , kTH 

R(J)=S(J+1) 

14 CONTI ilOF 

C ZKaUHTRUT OF KEBOUACK FILTFR 

7,K=0 

DO IS T al,JB 
7K»ZKtSC.I)*B(J) 

15 CIJNTTURE 
55COHLK = YK-ZK 

TF(5CURLK»GE*0) GO TO 25 

S(JB+l)s~t 
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r> f.iiHi isucounsTtOiCTEP h/\ta afire dee 

«m 'tn ih 

S(!-if i )s| 

TK (N. tf'TX) GO Tn 2,i 
TRfn.-ir. f £TUT+1 )) (,Yt I'P 2? 

FK =(SC»”iT.K-S( JBM. ) ) 

RK=L'' w wnn •> T GOAL 

» p ‘ 1r '' vv.vk of k as berimed farltar plus i after taps are fixed 
give i t i*: r i \ t i o j\i ,n ,T Yi fix taps gain 
TFfu. !;»:.•> ) r,*j to 22 

F‘ ! ItjT 1»|>JF * AGRA 

OH U <1 = 1 , JR 

R(.IJ = i} r >l)ArRU TAP GAIN VALUE Of FORWARD FILTER 

r(J)=fc n-F>mT[*f r ,n 

TYPE * » R ( t ) 

C.t'm'JUR 

DO lft 1 = 1 r.IH 

R C I J = UP.) ATR0 TAP GAIN VALUE OF FEEDBACK FILTER 
W(.U»B(JUFMI1LTI*SCJ) 

r ( J W T I i ! o E 

TRf SCJR+t ) .GT.O) GO TO 10 

TCK = -1.o . 

r,n T n Pli 
TCK = t 

TKCAfn .E0.1CK) GO TO 21 

D(T ) = '* # 

Mf)?=Ni)2ft 

WRITE C AG, t)W02,K 

TRfLn2.GT.100) GO TO 70 

on TO 28 

D(T}=' . ' 

GIVE SAME RANGE TO WRITE THE TAP GAIN VALVES IN 0/P FILE 
TFfK.LE. TRANCE) GO TO 27 
TFf K.GT.-7RANGE ) GO TO 27 
WRlTEny*85)CYK,ZK,SCURLK f SCJB+i) ,ECK)) 

WRITE ( 45>, 86 )(FF(.J) , J»1 , JF) , (BBC J) ,J=1,JS) 
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«’ OH'h '"slTW, UniaKrt UF F Vi 0 K BACK WARD TAPS 

P mm.mm.-v.m.m.-m.mm.'m'mm.mm.mmmmmmm-mm'mmmmmmmmmmm—m'mmm. 

Ab f 1 1 1 > : . A T f 2 X , 7 C E 1 3 . 4 ) ) 

77 i-.h r i !| :; 

io ru'M’i.O’E 

to i-notF'i ls,n >JF 1 

f »l K»x, 'NU OF ERROR T.RSS N0IRE= , »l7*5j{»f5*t) 

WHn-KOtb, J?)NlUS,ri02,K 

0 i FOP lATfSX, 'MOOBER OF ERROR AFTER SMCER AKE=% T 6 , 5X , ' NUMBER OF 

1 F.imjU IN l/P-O/R ARE= # , f4,*sx, 'VALVE OF K='I7) 

pf» R 4 [ = t,<)R 

C /PlT|.;f 1S,77HC lC(T*23-25+J) , IRCmWS+J) ,A(I*25-?5+J))»J=l»25) 

WHITE l lb, 37) CiiCl*iOO-100+J),.Itfl , 100 ) 

74 roRTiiOo: 

77 FOPOATf 1 OX, 100 (Al)) 

WRITE (4b,*)JEt 
W«TTFC4t>,*)NniS 
no os 1=1 , wr 

WKT'fF ( 4b, 97 ) (MC 1*100-100 + 0) ,0 = 1 ,100) 

OS ri.lNTTuHh 

WPTTKC45 ,09) 

Q9 FOP, iPTCtOX, 'ERROR IM ITERATION * ) 

no yl 1 = 1,2 

Wi/T’tT( 1S,RR)(E(I*10-t ) + J) ,J=1 ,10) 

Qt r 1 JOT TO' IE 

9 8 FORMAT ( 10( 1X,E12„4) ) 

STOP; FOR 



